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Abstract 

This communication presents the thermoeconomic optimization of an irreversible Stirling cryogenic refrigerator cycle. The 
external irreversibility is due to finite temperature difference between the working fluid and the external reservoirs while the 
internal irreversibility is due to the regenerative heat loss. The thermoeconomic function is defined as the cooling load divided 
by the total cost of the system plus the running cost. The thermoeconomic function is optimized with respect to the working 
fluid temperatures and the values for various parameters at the optimal operating condition are calculated. The effects of 
different operating parameters on the performance of the cycle have been studied. It is found that the effect of regenerative 
effectiveness is more than those of the other parameters on all the performance parameters of the cycle, for the same set of 
operating condition. 

© 2004 Elsevier Ltd and HR. All rights reserved. 
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1. Introduction 

The Stirling cycle is one of the important models of 
refrigeration systems for the production of very low 
temperature which has promoted the new developments in 
the design of this cycle. It is also desirable to have the 
corresponding results [1] for the coefficient of performance 
(COP) of the refrigeration systems. Leff and Teeters [2] 
have noted that the straight-forward Curzon-Ahlborn [1] 
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calculations will not work for a reversed Carnot cycle 
because there is no ‘Natural Maximum’ in these systems. 
Blanchard [3] has applied the Lagrangian method of 
undetermined multiplier to find out the COP of an 
endoreversible Carnot heat pump operated at minimum 
power input for a given heating load. 

In recent years the performance of the Carnot, Brayton, 
Stirling, Ericsson and other refrigeration cycles has been 
investigated by a number of researchers using the concept of 
finite time thermodynamics [4-9], the ecological criteria 
[10-17] and thermoeconomic approach [18-21] for 
different operating conditions. 

In this paper we will investigate the thermoeconomic 
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Nomenclature 

A 

area (m 2 ) 

C 

heat capacitance rate (kW K -1 ) 

Cf 

molar specific heat (kJ kmol -1 K -1 ) 

COP 

coefficient of performance 

P 

power (kW) 

Q 

heat transfer (kJ) 

R 

cooling load (kW) 

Ro 

gas constant (kJ kg -1 K -1 ) 

T 

temperature (K) 

t 

time (s) 

U 

overall heat transfer coefficient 


(kW m~ 2 K -1 ) 

1,2, 3, 4 

state points 

Greek 


s 

Effectiveness 

A 

Volume ratio 

Subscript 

a 

given in Eq. (20) 

b 

given in Eq. (22) 

c 

source-side 

f 

fluid 

h, H 

sink-side/heating 

L 

heat source/cooling 

max 

maximum 

m 

optimum 

R 

regenerator 


optimization along with a detail parametric study of an 
irreversible Stirling cryogenic refrigerator cycle for differ¬ 
ent operating conditions using the concept of finite time 
thermodynamics [1]. 


2. System description 

It is well known that the working substance of the 
Stirling cycle may be a gas, a magnetic material etc. For 
different working fluids, this cycle has different performance 
characteristics. When the working substance of the cycle is 
the prefect/ideal gas the cycle consists of two isothermal and 
two isochoric processes as shown on schematic and T-S 
diagrams in Fig. 1. This cycle approximates the expansion 
stroke of a real cycle as an isothermal process 1-2 with an 
irreversible isothermal heat addition at temperature T c from 
a heat source of finite heat capacity whose temperature 
varies from T Ll to T L2 . The heat addition to the working 
fluid from the regenerator is modeled as an isochoric process 
2-3. The compression stroke is modeled as an isothermal 
process 3-4 with an irreversible heat rejection at temperature 
T h to the heat sink of finite heat capacity whose temperature 



▲ 



cycle. 

varies from T m to T m . Finally, the heat rejection from the 
working fluid to the regenerator is modeled as an isochoric 
process 4-1, thereby completing the cycle. 

As mentioned earlier, the heat transfer processes 1-2 
and 3-4 in a real cycle must occur in finite time. This 
requires that these heat processes must be proceeded 
through a finite temperature difference and therefore, 
defined as being externally irreversible. There is also some 
heat loss through the regenerator, as an ideal regeneration 
requires an infinite regeneration time or area, which is not 
the case in practice. 

3. Thermodynamic analysis 

Let Q c be the amount of heat absorbed from the source at 
temperature T c and Q h be the amount of heat released to the 
sink at temperature T h during two isothermal processes, 
then: 

fih = T h AS = C H (Xm ~ Tm)tn = tf H A H (LMTDj H f H (1) 
Qc = T C AS = C L (T Ll - T L2 )t L = t/ L A L (LMTD) L / L (2) 
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where A5(= nR 0 In A) is the entropy change, A is the 
volume ratio, n is the number of mole for the working fluid, 
R 0 is the universal gas constant, C L , C H and t L , t H are the 
heat capacitance rates of source/sink reservoirs and heat 
addition/rejection times, s H and s L are the effectiveness of 
the heat exchangers for hot- and cold-side, respectively, 
f/ H A H and U L A L are the overall heat transfer coefficient-area 
products and (LMTD) H and (LMTD) l are the log mean 
temperature difference on sink- and source-side, respect¬ 
ively, defined as: 

(LMTD)h = [<T h - r H1 ) - (T h - T m )]/ln[(T h 

- T m )/(T h - T m )\ (3) 

(LMTD) l = [<T U - T e ) ~ (T L2 - r c )]/ln[(r L1 


absorbed from the source are: 


e H = G h -Ag R (ii) 

Ql = Qc - Ag R (12) 


For a Stirling refrigeration cycle, the irreversibilities 
associated with the heat transfer during the compression and 
expansion of the working fluid are a significant source of 
thermodynamics inefficiencies of the system. Thus we have: 


Qh 

Th 



(13) 


where f H and T L are, respectively, the average temperatures 
of the sink- and source-fluid at constant pressure and 5 gen is 
the entropy generation rate of the cycle. Using Eqs. (1), (2), 
(11) and (12), we have: 


- T c )/(T L 2 ~ TJl 


(4) Q h /T h - Q c /T c = 0 


Using Eqs. (1) and (2) and Eqs. (3) and (4), we can get: W = Q n — Q L = Q h — Q c 


(14) 

(15) 


Qh = ^H s H(Th — T m )t u (5) 

Qc — C L s L (T Ll — T c )t L (6) 

where s H = 1 — e -(c/A)H/CH and s L = 1 — e _(t/A)L/CL are the 
effectiveness of the hot- and cold-side heat exchangers, 
respectively. 

As these cycles, in general, do not possess the condition 
of perfect regeneration, it is assumed reasonably that the 
regenerative heat loss per cycle is proportional to the 
temperature difference of the two isothermal processes as 
assumed by earlier workers [5-9]. 

A£ r = nc f (l ~ s R )(T h - T c ) (7) 

where c f is the molar specific heat capacity of the working 
fluid and s R is the effectiveness of the regenerator, defined 
as: 

^ Qregen ,actual QaX QlY = T b -T x 

(2regen,ideal Qa I QlS -^h T c 


1 Y — * r 2 Vr 

= —-- - -( 8 ) 

T h -T c 1+A r 

where N R = (UA) R /Cf is the number of transfer unit of the 
regenerator and C f is the heat capacitance rate of the 
working fluid. 

Owing to the influence of irreversibility of finite heat 
transfer, the regenerative time should be finite as compared 
to that of the two isothermal processes [5-9,17]: 

Ir = h + h — 2tr(r h — T c ) (9) 

where a is the proportionality constant, which is indepen¬ 
dent of the temperatures on the hot- and cold-side but 
depends on the property of the regenerative material. Thus, 
the total cycle time £ cycle will be: 

t C yc\e = (h-r + A + A) (10) 

When the irreversibilities mentioned above are taken 
into account, the net amounts of heat released to the sink and 


It is worthwhile to note that for an irreversible Stirling 
cycle Eq. (14) is true, although the irreversibilities 
associated with the heat transfer are considered. 

Since, the power input, cooling load and COP are the 
important parameters of a refrigerator cycle. Using Eqs. 
(11)—(15), we have: 

D _ Ql _ Qc ~ AQr 

AI — - — - 

Cycle % + 


[1 -fljfct- 1)] 


[ 




+ 


1 


p = 


(xy T m ) k L (T Ll y) 
W (Q h - Q c ) 


+ b i (pc 1) 


] 


(16) 


t. 


cycle 


? H + ? L + f R 

(x- 1 ) 


[ 




+ 


1 


+ bi(x — 1) 


COP = 


(xy C H] ) k L (T L{ y) 

Rl Qc~^Qr [l-a^x-l)] 


} 


(17) 


Qh Qc 


(x-l) 


(18) 


where k±i — s^Cn, k R — s R C R , x — T^/T c , y — T c , b\ — 
2al(nRQ In A) and = c f ( 1 — s r )/Rq In A. 

The objective function of thermoeconomic optimization 
proposed by earlier workers [18-21] is given by: 


Ql 

Q + c e 


(19) 


where Q and C e refer to annual investment and energy 
consumption costs, respectively. The investment cost was 
considered the costs of the main system components that are 
the heat exchangers and the compression and expansion 
devices together. The investment cost of the heat exchangers 
is assumed to be proportional to the total heat transfer area 
[18-21]. On the other hand, the investment cost due to the 
compression and expansion devices is assumed to be 
proportional to their compression/expansion capacities or 
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the required power input. Thus the annual investment cost of 
the system can be given by: 


C\ — a(A n + A l + A R ) + b 2 P 


— a(A n + A l + A r ) + b 2 (Qh GcVCycle (20) 

where the proportionality constant for the investment cost of 
the heat exchanger, a is equal to the capital recovery factor 
times investment cost per unit heat exchanger area and 
the proportionality constant for the investment cost for the 
compression and expansion devices, b 2 is equal to the 
capital recovery factor times investment cost per unit power 
input. The initial investment cost is converted to equivalent 
yearly payment using capital recovery factor [18-21]. The 
annual energy consumption cost is proportional to the power 
input. 

Q = ^3^ = ^(Gh — GcVCycle (21) 


where the coefficient, b 3 is equal to the equivalent annual 
operation hours corresponding to power input times price 
per unit energy [18-21]. Substituting Eqs. (20) and (21) into 
Eq. (19), we have: 


_ Ql _ 

a (An +^L +^r) + KQh ~ GcVCycle 


( 22 ) 


where b = b 2 + b 3 . Thus from Eqs. (23), (16) and (17), we 
have: 


bF = 


[1 - «x(jc - 1)] 


k x k 


[ 


X 


+ 


1 


^h(T ) 7 ^Ht) k l (t u y) 


+ bi(x— 1) j + (x - 1) 


(23) 


where k { = alb and 

Q £ r Qi 


ki = 


( 


f/ R (l - C R ) U 


H 


ln(l e H ) —— ln(l e L ) 

c 7 t 


) 


It can be seen from Eq. (23) that bF is a function of two 
variables x and y (as other parameters are constant for a 
typical set of operating conditions). Thus optimizing bF 
with respect to ‘y’ viz. dbF/dy = 0, yields: 

x (T L i — y) = y/ku/k L (xy — T H1 ) (24) 

Substituting Eq. (24) into Eqs. (23), (16) and (17), we have: 
[1 - fll (jc-l)] 


bF (x ) = 


xk< 


(25) 


R Ux) ~ 


(xT Ll r H1 ) 

1 — ai(x — 1 ) 


+ b 4 (x - 1) 


xk, 


(26) 


P(x) - 


(r H i ^r L1 ) 

(X- 1) 


+ by{x — 1) 


xk, 


(27) 


(Thi xTli) 


+ bi(x— 1) 


where k 3 = ^Jk u /k L , k A = (1 + k 3 )(k H l + k L l k 3 ! ), k 5 = 
k { k 2 k 4 , and b 4 = 1 + k { k 2 b { . 


4. Discussion of results 

In order to have the numerical appreciation of the results 
for the Stirling refrigerator cycle, we continue to investigate 
the effects of the effectiveness of different heat exchangers 
(e H , s L and s R ), the economic parameter (k { ) and heat 
capacitance rates (C H and C L ). The effects of each one of 
these parameters are examined while the rest of the 
parameters are kept constant as (e H = s L = £ r = 0.80, 
T m = 290 K, T L1 = 200 K, k { = 0.50, A = 2.0, 
C U = C L = 1.0 kW/K, b { = 0.01, U H = U L = U R = 2.0 
kWnr), the results obtained are as follows. 

4.1. Effects of cycle temperature ratio (x) 

The variation of the objective function, COP, power 
input and cooling load with respect to the cycle temperature 
ratio ( x = T h /T c ) for a typical set of operating parameters is 
shown in Fig. 2. It is seen from Fig. 2 that the objective 
function and the cooling load first increase and then 
decrease while the COP monotonically decreases whereas 
the power input monotonically increases as the cycle 
temperature ratio increases. These properties may be 
directly expounded by Eqs. (25)-(27) and (18), because 
the objective function and cooling load are not mono tonic 
functions of x while the power input and COP are 
respectively, the monotonically increasing and decreasing 
functions of x. It can also be clearly seen from Fig. 2 that 
both the objective function and the cooling load attain their 
maxima but at different value of x and there exists the 



Fig. 2. Effects of cycle temperature ratio (x = T h /T c ) on the various 
performance parameters of an irreversible Stirling refrigerator 
cycle. 
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Fig. 3. (a)-(d) Effects of various effectiveness on the maximum objective function and the corresponding COP, cooling load and power input of 
an irreversible Stirling refrigeration cycle, (a) Maximum objective function vs. effectiveness, (b) Corresponding COP vs. effectiveness, (c) 
Corresponding cooling load vs. effectiveness, (d) Corresponding power input vs. effectiveness. 


following relation: 

(.X 0 pt) bF — C^opt)i? L (28) 


where (x opt ) hF and (x opt )^ L represent the two different 
optimal values of x, the former one corresponds to the 
point of the maximum objective function while the later one 
corresponds to the point of the maximum cooling load, 
respectively. It is seen from Eqs. (25) and (26) that both the 
performance parameters, i.e. bF and R L are the functions of 
a single variable, v, for a typical set of operating condition. 
Thus, maximizing bF and R L with respect to x, yields: 




L 


B ± yJ(B 2 -AC) 

A 

B\ — -\J(B 2 — A t Ci) 
^1 


(29) 

(30) 


where A = [(b^T^i + B = b^T^Tm, C = 

{^4^hi ~~ ^5(1 + a \))Tmi ^1 — [(^i^li +a 1 & 4 ]r L1 , B { = 
biT hl T m and C\ = {b { T m — K 4 ( 1 + ai))T m . 


It can be seen from Eqs. (29) and (30) that both the 
parameters, i.e. (v opt )^ F and (v opt )^ L , have two roots for a 
given set of operating parameters, but the roots containing 
the minus sign do not give the useful results. So we 
calculated the values of (v opt )^ and (v opt )^ L only with the 
plus sign. Substituting the values of (v opt )^, F into Eqs. (25)- 
(27) and (18) we can calculate the maximum objective 
function and the corresponding cooling load, power input 
and COP, while the maximum cooling load and the 
corresponding power input and COP can be calculated by 
substituting (v opt ) P into Eqs. (18) and (27). On the other 
hand, the optimal values of y, for both the cases, i.e. Cy 0 pt)z?F 
and Cy 0 pt)# L can be calculated by substituting the values 
(a opt)bF an d (^ 0 pt)i? L separately, into Eq. (24) for a typical set 
of operating condition. 

4.2. Effect of effectiveness 

Fig. 3(a)-(d) show the effects of the effectiveness on the 
maximum objective function, the corresponding COP, 
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cooling load and power input of an irreversible Stirling 
cryogenic refrigeration cycle. It is seen from these figures 
that the maximum objective function and the corresponding 
COP first increase and then decrease while the correspond¬ 
ing cooling load and power input increase as the effective¬ 
ness on any heat exchanger increases. Their physical 
meaning may be explained as follows. For the typical set 
of operating parameters given above, it implies the fact that 
the larger the effectiveness is, the larger the heat transfer 
areas are required and the smaller the irreversibility 
associated with the cycle is. This results in the larger values 
of the corresponding cooling load and power input. In 
general, the corresponding cooling load is not a linear 
function of the corresponding power input, so that there may 
be a maximum for the COP. On the other hand, when the 
heat transfer areas are increased, the corresponding cooling 
load and the cost of the system increase. In general, the 
corresponding cooling load is not a linear function of the 
cost of the system, so that there may be a maximum for 
the objective function. As —► 1, A, —*■ oo (I = H, L, R) so 
the size of the system and hence, the cost becomes very 
large. This results a sharp increase in the corresponding 
cooling load and power input while a sharp decrease in the 
objective function and the corresponding COP as can be 
seen from Figures 3(a-d). Again, as the effect of the internal 
irreversibility is more than that of the external one and the 
regenerative-side effectiveness belongs to the internal one 
unlike the hot- and cold-side effectiveness, which belongs to 
the external one. Hence, the effects of the regenerative 
effectiveness are more than those of the hot- and cold-side 
effectiveness on all the performance parameters for the same 
set of operating condition. Moreover, the effects of the hot- 
and cold-side effectiveness are almost the same for all 
performance parameters, so these two curves overlap, as can 
be seen from these figures. 

4.3. Effect ofk { 

The effects of economic parameter (k { ) on the maximum 
objective function, the corresponding cooling load, power 
input and COP are shown in Fig. 4(a) and (b). It is seen from 
these figures that the maximum objective function and the 
corresponding COP decrease while the cooling load and the 
power input increase as the economic parameter (kf) 
increases. Since, the cost of the system increases by 
increasing the economic parameter so the maximum 
objective function decreases with increasing the economic 
parameters. Also the optimal value of cycle temperature 
ratio (x opt ) bF increases with increasing k x , resulting in an 
increase in the corresponding cooling load and power input 
but a decrease in the objective function. Since, the 
corresponding power input in this region increases more 
than that of the corresponding cooling load, hence, as a 
result the COP decreases, by increasing the economic 
parameter, which can be seen clearly from these figures. 




Fig. 4. (a)-(b) Effects of the economic parameter (k\) on the 
maximum objective function and the corresponding COP, cooling 
load and power input of an irreversible Stirling refrigeration cycle, 
(a) Maximum objective function and optimal COP vs. economic 
parameter, (b) Corresponding cooling load and power input vs. 
economic parameter. 

4.4. Effect of heat capacitance rates 

The effects of source/sink-side heat capacitance rate on 
the maximum objective function, the corresponding COP, 
cooling load and power input are shown in Fig. 5(a) and (b). 
It is seen from these figures that the maximum objective 
function as well as the corresponding COP decrease, as the 
heat capacitance rate on source/sink-side reservoir 
increases. On the other hand, the corresponding cooling 
load or power input is found to be a monotonically 
increasing function of the heat capacitance rate on the 
source/sink-side. These results can be explained as follows. 

For the typical set of operating parameters given above, 
it implies the fact that the larger the heat capacitance rates 
are, the larger the heat transfer areas are required. This 
results in the larger values of the power input, cooling load 
and the cost of the system. As a result the maximum 
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(b) C H or C L (kW/K) 

Fig. 5. (a) and (b) Effects of the source/sink-side heat capacitance 
rate (C H or C L ) on the maximum objective function and the 
corresponding COP, cooling load and power input of an irreversible 
Stirling refrigeration cycle, (a) Maximum objective function and 
optimal COP vs. heat capacitance rate, (b) Corresponding cooling 
load and power input vs. heat capacitance rate. 


objective function and the corresponding COP decrease 
while the cooling load and power input increase, as the heat 
capacitance rates of the external fluid increases. 


5. Conclusions 

A more realistic Stirling refrigerator cycle model 
including external and internal irreversibilities for the finite 
heat capacities of external reservoirs has been studied in 
detail. The external irreversibilities are due to finite 
temperature difference between the cycle and the external 
reservoirs, while the internal irreversibility is due to the heat 
loss inside the regenerator. The thermoeconomic function, 
which is defined as the cooling load per unit cost of the 


system (viz. running, investment etc.) has been adopted as 
an objective function for maximization. The objective 
function is maximized with respect to the cycle tempera¬ 
tures and the corresponding cooling load, power input and 
COP are evaluated for the different operating conditions. 
The objective function is found to be a decreasing function 
of the economic parameter and heat capacitance rate while 
there are optimal values of the cycle temperature ratio (x) 
and different effectiveness, at which the objective function 
attains its maximum value for a typical set of operating 
condition. It is also found that effects of the regenerative- 
side effectiveness is more than those of the other side 
effectiveness, which can be explained in terms of internal 
and external irreversibilities. On the other hand, the effects 
of the source- and sink-side parameters are found to be 
similar, so the two curves overlap, as can be seen from 
graphs. Thus, the present cycle model gives some optimal 
criteria, which will be useful to understand and optimize the 
design performance of a real cycle from the point of view of 
thermodynamics as well as from the point of view of 
economics. The results obtained here are also applicable to 
the Ericsson cycle, in which the only difference is that the 
two isochoric processes will be replaced by the two isobaric 
processes and the volume ratio will be replaced by pressure 
ratio and hence, the specific heat at constant volume will be 
replaced by the specific heat at constant pressure. 
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